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Abstract. Elementary complex analysis and Hilbert space methods show 
that a union of at most n arcs on the circle is uniquely determined by the 
nth Fourier partial sum of its characteristic function. The endpoints of the 
arcs can be recovered from the coefficients appearing in the partial sum by 
solving two polynomial equations. 



We let T = {z e C : \z\ = 1} and D = {z G C : \z\ < 1}, and for any subset E of 
T and integer k we write 



for the fcth Fourier coefficient of the characteristic function xe of E. As bounded 
functions with the same sequence of Fourier coefficients agree almost everywhere, 
any subset E of T is determined up to a set of measure zero by the sequence E(k). 
If E is known to have additional structure, the entire sequence may not be needed 
to recover E. Our present subject is a simple yet nontrivial illustration of this 
principle. 

An arc is by definition a closed, connected, proper and nonempty subset of T. 
We declare T along with the empty set to be a "union of arcs." 

Theorem 1. If n is a nonnegative integer and E\ and Ei are unions of at most n 
arcs satisfying 



then Ei = Ei . 

Thus a set E that is known to be a union of at most n arcs can be recovered 
completely from the nth Fourier partial sum of xe, regardless of any quantitative 
sense in which this partial sum fails to approximate xe- This stands in slight 
contrast to the well-known defects of Fourier partial sum approximation of functions 
with jump discontinuities, such as the Gibbs phenomenon (see e.g. [H Chapter 17]). 
Significantly, the property of the Fourier basis expressed by Theorem[T]is not shared 
by other orthonormal systems of functions on T (see $3]). 

Our proof of Theorem [T] exploits a connection between unions of arcs and certain 
rational functions — the Blaschke products, whose properties we recall in $1] Each 
Blaschke product has a nonnegative integer order. In [21 we construct an injection 
E t->- bs from the set of finite unions of arcs to the set of Blaschke products with 
the property that if E is a union of at most n arcs, then oe has order at most n. 
This map has the property that if E\ and Ei satisfy ([1]), then oex and &e 2 have 
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the same nth order Taylor polynomial at 0. To prove Theorem Q] it then suffices to 
note, as we do in |j3l how a Blaschke product of order at most n is determined by 
its nth order Taylor polynomial. 

With Theorem [1] in hand, one may ask how to recover E from a partial list of 
Fourier coefficients in an explicit fashion. This is the subject of SJH where we present 
an algorithm for testing whether or not a given tuple of complex numbers takes the 
form (I?(fc))£_Q for a union E of at most n arcs, and for finding the endpoints of 
these arcs in terms of the Fourier coefficients in this case. 

Perhaps because of its elementary nature, we have not found Theorem[T]explicitly 
stated in the literature, although it is known, and the literature abounds with more 
general theorems on the reconstruction of a function from partial knowledge of its 
Fourier transform. In [6] it is shown that a function on T that is piecewise constant 
on a partition of T into m connected pieces may be recovered from its mth Fourier 
partial sum. Note that Theorem [T] concludes slightly more from a much stronger 
hypothesis. 

The argument we use is known to specialists. The basic idea is to apply a 
conformal map into the disc and then the classical Caratheodory-Fejer theorem pQ. 
This is by no means the only approach to Theorem [TJ It should be contrasted with 
what one may get by viewing ([T]) as a system of polynomial equations and solving 
it directly with algebra. 

We are indebted to Donald Sarason for many valuable discussions, and to Mihalis 
Kolountzakis for drawing our attention to [6]. 

1. Blaschke products 
Definition. A (finite) Blaschke product is a function of the form 

n 

( 2) >(^nS| 

j= i j 

for some nonnegative integer n, some A 6 T, and some a\,...,a n £ D. The 
nonnegative integer n is called the order of the Blaschke product. 

If n = we interpret the empty product as 1 . The domain of a Blaschke product 
is either T, D, or the closure D of D, depending on context. A Blaschke product is 
evidently a rational function that maps T to itself and has no poles in D (it suffices 
to check the case n = 1). It is well known that these properties characterize the 
Blaschke products. 

Proposition 1. If a rational function r maps T to itself and has no poles in D, 
then it is a Blaschke product of order equal to the number n of zeros of r in ID, 
counted according to multiplicity. 

Proof. We induct on n. If n = 0, then r = q^ 1 for some polynomial q; write 
l( z ) = SfcLo 1 kzk wrtn 9™ ^ 0. As q(T) C T we have 

m T-^m — m — k 

q (z)- 1 =W)= Wr 1 ) = E^~ fe = ,T — - * e T, 

so this holds for all nonzero z € B. As q has no zeros in D, the extreme right hand 
side has no pole at 0; thus m = and q is constant as desired. 

If r has n + 1 zeros in D, choose one, a, and note that r(z) ■ ( r5§z ) -1 has n zeros 
in D and maps T to itself. □ 
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Definition. If b is a Blaschke product, we let Ub = {z G T : hnz > 0}. 
If the zeros of a Blaschke product are oi, . . . , a n , we calculate from ^ 



zb'jz) 

K z ) j-^\ z ^ a 3 



El — a,- 



so the argument of b{e lt ) is strictly increasing in i. The argument principle implies 
that b(e tt ) travels n times counterclockwise around T as t runs from to 2ir. 

Corollary 1. A Blaschke product b has order n if and only ifUb is a disjoint union 
of n arcs. 

This is the main reason we include T as a "union of arcs." 

2. Blaschke products from unions of arcs 
Let 5 = {zeC:0<2Rez<l} and let <f> denote the function 

= exp(2m(z - 1/4)) - 1 
<PV "' exp(27ri(z - 1/4)) + 1' 

It is easy to show (see e.g. [U §111.3]) that 4> maps S bijectively onto D \ {±1}, 
that <f> restricts to an analytic bijection of the interior of S with D, that </> maps 
the right boundary line of S onto {z g T : Imz > 0}, and that <f> maps the left 
boundary line of S onto {z G T : Imz < 0}. 

Proposition 2. If E is a disjoint union of n > arcs and He is given by 

1 oo 

(3) h E (z) = -E{0) + ^2E(k)z k , zeO, 

fe=i 

then He is an analytic map o/D into S, and the function D — > D given by 

b E = 4>°h E 

extends uniquely to a Blaschke product D-)B o/ order n satisfying Ub E = E. 

Using the formulas for <f> and h E one can show without much work that b E is 
a rational function; the work in proving Proposition [5] is to establish that b E has 
the mapping properties of Proposition [I] and hence is a Blaschke product, and to 
prove that Ub E = E. 

To motivate the argument, let us work nonrigorously for a moment. Formally 
we have the series expansion 

(4) X e(z) = ^%)z fe , zeT, 

feez 

and formal manipulation of the series ([3]) with z G T then shows that 



Xe{z) = h E (z) + h E {z) = 2 Re h E (z), z G T. 

As \E is {0, 1} valued on T, the maximum principle for harmonic functions then 
implies that h E maps U into S, so b E ~ cf) o h E maps D into D and sends the circle 
to itself. By Proposition Q] it follows that b E is a Blaschke product; the equality 
Ub E — E comes from the mapping properties of <fi on the boundary of S. 

What makes this argument nonrigorous is that the series (j4]) does not converge 
for all z G T, and to equate xe with 2 Re h E is to ignore the distinction between 
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a discontinuous real valued function on T and a harmonic function on D. To fill 
in these gaps, we need to use the actual connection between 2 Re h E and xe — the 
former is the Poisson integral of the latter. 

Proof. It is easily checked that © does define an analytic function on D, e.g. 
because 2fc=i l-^(^)! 2 is convergent. One can then verify the identity 

2h E (z) = — / — — -XE(e ls ) ds, z G D. 

^ 7r Jo 1 ~ ze 

(Fix z, expand 1 _ z 1 e _ is as a power series in z and interchange the sum and the 
integral.) Taking real parts it follows that for any r G [0, 1) and any t 

1 C 2n 

(5) 2Reh E (re lt ) = — / P r (t - s) XE (e ls ) ds, 

where 

P r (t) = Re 



1 



1 — re 

is the Poisson kernel. It is elementary (see e.g. [3J §X.2]) that for r G [0, 1) the 
function P r is nonnegative and satisfies ^- J Q 27r P r (8) d6 — 1; thus ([5]) implies that 
2Re h E (z) G [0, 1] for all z G D, and /i_e maps D into 5. 

As r increases to 1, the P r converge uniformly to the zero function on the com- 
plement of any neighborhood of (see e.g. [5J §X.2]). From we conclude 

(6) \im2Re h E (rz) = xe(z) 

rfl 

at any z G T at which XE is continuous. We conclude that for any such z the limit 
lim r -|-i((/> o h E )(rz) exists and is in T. 

We claim that <j) o h E is a rational function. In the case n = this is clear. 
Otherwise, from the definition of <fi it suffices to show that exp(2i:ih E ) is a rational 
function, and for this it suffices to treat the case n = 1. In this case there are 
real numbers a < b with b — a < 2ir satisfying E = {e lt : t G [a, b}}, and i?(fc) = 
exp ( - 1 kb 2~°*v ( - 1 fc a I for all fc > 0. Let log denote the analytic logarithm defined 
on C \ {z G C : z < 0} that is real on the positive real axis and recall that 
log(I — z) = — J2kLi X f° r au 2 G B. A comparison of power series shows 

h E (z) = ^ + -L (log(l - e- 4fc z) - log(l - e—z)) , z G D, 

so exp(27ri/i£;) = exp(i^^) \z~^%sr is rational. 

At this point we know that b E = <p ° is a rational function mapping D into 
itself. From (|6]) we deduce that b E maps T into itself, so b E is a Blaschke product 
by Proposition [T] The equality Ub E — E then follows from ([6|). The order of b E is 
n by Corollary [T] □ 

If E\ and are two unions of arcs related by (jTj), it is clear from the definition 
that h El and h E2 have the same nth order Taylor polynomial at 0. As <p is analytic 
at 0, the same is true of b El and b E2 . 

Corollary 2. If n > and E\ and Ei are each unions of at most n arcs satisfying 

(7) E\(k)=E2(k), 0<k<n, 
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then there are Blaschke products b\ and hi, each of order at most n, satisfying 
Ej = Ubj for j — 1,2 and 

(8) b 1 {k)=b 2 {k), 0<k<n. 

3. Blaschke products from Toeplitz matrices 

Fix a positive integer n for the remainder of this section. Our goal is to show 
that Blaschke products b\ and 62 having order at most n and satisfying ([SJ must 
be equal. Let L 2 denote the space of square-integrable functions T — 5- C, with inner 
product 

</,<?) = ^ / " f(e U )Wt)dt, f,g e L 2 . 
27T Jo 

(We identify two functions if they agree almost everywhere.) 

For < k < n we let £ fe denote the function T — ► C given by z i-> z k . It is 
immediate that {( k : < k < n} is an orthonormal subset of L 2 . We denote its 
span, the space of analytic polynomials of degree at most n, by P; we let ir : L 2 — > P 
denote the orthogonal projection. 

Definition. If f : T — > C is bounded, Tf : P — » P denotes the linear map given by 

TfS, - tt(/0, e e P 

Here /£ is t/ie pointwise product of f and t;. 

If we let ||T/|| denote the norm of T/ regarded as a linear operator on P and 
write ll/Hoo = sup 2gT it is clear that 

\\Tf\\ < ll/lloo 

for any bounded /. It is also clear that for any such / 

(T f t k X j ) = fU-k), 0<j,k<n, 

so the matrix of Tf with respect to the orthonormal basis {( k : < k < n} is 
constant along its diagonals (it is a Toeplitz matrix). 

If / is a Blaschke product, then / is analytic on D, so the matrix of Tf is lower 
triangular with first column (f(k))^ =0 . Our hypothesis (|8]) is thus that = Xb 2 , 
and to deduce that 61 = 62 it suffices to show how to recover a Blaschke product b 
of order at most n from the operator Tb it induces on P. 

Lemma 1. If b is a Blaschke product of order at most n, then ||Tfc|| = 1, and for 
any nonzero r £ P satisfying ||Tf,r|| = ||r|| one has T^r = br. 

This proof is a special case of the proof of [51 Proposition 5.1]. 

Proof. There are nonzero polynomials p and q, each of degree at most n, satisfying 
b = p/q. Clearly Ti,q = p, and as b maps T to itself, we have \p(z)\ = \q(z)\ for all 
z e T, so ||p|| = ||g||. We deduce that ||T 6 g|| = ||<?|| and thus ||T b || > 1; since also 
\\T b \\ < H&lloo = 1, we conclude ||T 6 || = 1. 
If r G P satisfies \\Ti,r\\ = \\r\\ we have 

||r|| a = ||T 6 r]] 2 = ||7r(6r)|| 2 < ||6r]] 2 = [*" \b(e*)\ 2 \r{e u )\ 2 dt = ||r|| 2 , 

Jo 

from which ||7r(6r)|| = ||6r|| and thus 7r(6r) = br as desired. □ 
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Remark 1. The argument of Lemma\]\can be modified to show that if f is bounded 
and analytic on D and ||/||oo = 1, then < 1 with equality if and only if f is a 
Blaschke product of order at most n. With more work, one can prove the rest of the 
classical Caratheodory-Fejer theorem: that every lower triangular (n + 1) x (n + 1) 
Toeplitz M satisfying \\M\\ = 1 is of the form Tf for such an f . 

We can now prove Theorem [TJ 

Proof of Theorem[J\ By Corollary [2] there are Blaschke products b\ and bi of order 
at most n satisfying U bj = Ej for j = 1,2 and bi(k) = ^(fc) for < k < n. This 
second fact implies that T bl = T b2 . By Lemma [T] there is nonzero q s P satisfying 
\\T bl q\\ = \\T b2 q\\ = ||<?|| and 

T bl q T b2 q 

Oi — — = 02, 

so ^ = U bl = U b2 = E 2 . □ 

As the Fourier coefficients of a bounded function are coefficients with respect to 
an orthonormal basis of the Hilbert space L 2 , one might wonder if Theorem [T] is a 
special case of a simpler result about arbitrary orthonormal bases of L 2 . This is not 
the case. There are, for example, orthonormal bases B for L 2 with the property 
that for every finite subset FCB, there is an arc A with the property that every 
element of F is constant on A. (The basis (e 27rit M> f{t))feH, where H is the 
Haar basis of L 2 [0, 1] constructed in [3l §111.1], has this property.) In this situation, 
if E C A and E' C. A are any two unions of arcs with the same total measure, 
one will have (xE,f) = (XE'if) f° r au / S F: any finite collection of coefficients 
with respect to B must fail to distiguish infinitely many unions of n arcs from one 
another. 

4. An algorithm 

Let J- denote the map sending a union of at most n arcs E to the tuple (E(k))^ =0 
in C" + . Suppose c = (ck)k=a ^ s gi yen : an d we desire to know whether or not c 
in the range of J- ' . The arguments of the previous sections give us the following 
procedure. (We use the orthonormal basis of to identify linear operators on P 
with (n + 1) x (n + 1) matrices.) 

(1) Calculate the nth Taylor polynomial at for 4>{^ + Sfe=i c kZ k ), and make 
its coefficients the first column of a lower-triangular Toeplitz matrix M. 

(2) Evaluate ||M||. 

If 1 1 A/ 1 1 7^ 1, then c is not in the range of T . 

(3) Otherwise ||M|| = 1 and by the Caratheodory-Fejer theorem (see Remark[T]) 
there is a unique Blaschke product / of order at most n satisfying M = Tf. 
Find F — Uf (e.g. by solving f{z) = ±1 to get the endpoints of the arcs) 
and calculate the coefficients of the nth order Taylor polynomial at for 
b F . 

If these coefficients are the first column of M then bp = f and c = J-"(F); 
otherwise c is not in the range of T . 

Remark 2. The third step of the algorithm is necessary as the map E M> bp from 
unions of n arcs to Blaschke products of order n is not surjective. One can check, 
for example, that of the Blaschke products bt(z) — fcrBl for real \t\ < 1, all of 
which satisfy U bt — U b „, only bo is in the range of E h-> bp. 
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If we know in advance that c = J-(E) is in the range of J 7 , this algorithm can 
recover E from c in a somewhat explicit fashion. The matrix M constructed from c 
is Tft E ; Lemma [1] implies that if we choose a nonzero q 6 P satisfying ||-Mg|| = 
we will have 6e = If q is chosen so as to have minimal degree, the polynomials 
Mq and q will have no nontrivial common factors. In this case the degree of q is 
the order of bg, and the endpoints of the arcs of E — the solutions to 6e(z) = 1 and 
i>E{z) — — 1 — are the roots of the polynomials Mq — q and Mq + q. A computer has 
no difficulty carrying out this procedure to find the arcs of E to any given precision 
from the tuple c = F(E) . 

As this algorithm involves solving polynomial equations, we cannot expect sym- 
bolic formulas for these endpoints of the arcs of E in terms of the Fourier coefficients 
E(k). Formulas for the polynomials Mq ± q, however, can be obtained with some 
effort. The entries of M are polynomials in exp(27rzi?(0)), E(l), E{n) with 
complex coefficients. As M has norm 1, a vector q will satisfy = ||g|| if and 

only if q is an eigenvector for the self-adjoint matrix M*M corresponding to the 
eigenvalue 1; we can find such a q by using Gaussian elimination, for example. As 
the entries of M* M are polynomials in the entries of M and their complex conju- 
gates, the coefficients of q and Mq ± q will be rational functions in exp(27ri£'(0)), 
E{1), E{n) and their complex conjugates. Cases may arise in computing 
Mq ± q symbolically: in row reducing the symbolic matrix M*M — J, one needs to 
know whether or not certain functions of the matrix entries are zero — but explicit 
formulas can be obtained in every case. 

We give one example. Suppose that E is a union of at most two arcs, with E(0), 
E(l), and E(2) given. Write E = exp(2iriE(0)) and E k = -2nikE(k) for k = 1, 2. 
Carrying out the above procedure, one finds that if both E\ and the denominator 

of _ _ 

E<iE\ ~h 2,Ei — E^Ei — 1E\E{) 

E\Eq + E 2 E Q -E 2 + El ' 
are nonzero, then the starting points of the arcs of E are the solutions z of the 
equation 

a ( E-i + {l- E )a \ 

Z ~ aZ+ { E lEy J=°- 
The endpoints of the arcs of E are given by a similar formula. 
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